Passive Sliders on Fluctuating Surfaces: Strong-Clustering States 
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We study the clustering properties of particles sliding downwards on a fluctuating surface evolving 
through the Kardar-Parisi-Zhang equation, a problem equivalent to passive scalars driven by a 
Burgers fluid. Monte Carlo simulations on a discrete version of the problem in one dimension reveal 
that particles cluster very strongly: the two point density correlation function scales with the system 
size with a scaling function which diverges at small argument. Analytic results are obtained for the 
Sinai problem of random walkers in a quenched random landscape. This equilibrium system too has 
a singular scaling function which agrees remarkably with that for advected particles. 
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Interesting correlations develop in a system of particles 
which are driven by a fluctuating field but do not react 
back on it. This problem of semiautonomously coupled 
systems is of general interest and importance. In the 
context of fluid dynamics, it comes up as the passive 
scalar problem, where a density field (e.g. dye particles) 
or a temperature field is advected by a stirred fluid 0,0- 
The interplay of fluid advection and random thermal 
noise can lead to two distinct behaviours of the passive 
scalar field — it may either spread out Q, [j| , or show 
a tendency to collapse and form clustered aggregates. 
The clustered state can arise when either the passive 
scalar flow p| or the fluid itself 5| is compressible. 
Clustering also arises in models of thin film growth in 
binary systems where domain walls between the two 
different species can act like passive particles while 
surface growth is the driving field 0, Q • 

In this paper, we study the clustered state that 
occurs in a simple but nontrivial model, namely non- 
interacting, passive particles sliding downwards on a 
fluctuating one-dimensional interface described by the 
Kardar-Parisi-Zhang (KPZ) equation. This problem can 
be mapped onto that of particles advected by a noisy 
Burgers fluid. Some aspects of this problem have been 
earlier studied by Drossel and Kardar 0,0- However, 
the questions we study here are different and our results 
indicate that the system reaches a new type of steady 
state, namely the strong-clustering state (SCS). The 
signature of the SCS is the divergence at origin of 
the two-point density-density correlation function as 
function of the separation scaled by the system size. 
This divergence provides a measure of clustering in the 
system, and distinguishes the SCS from normal phase 
separated states || and also from fluctuation-dominated 
phase separated states 0, ^3 where a similar scaling 



function approaches a constant value, though in a 
singular fashion, at small argument. 

Our numerical results show that the SCS occurs over a 
wide range of parameter values in the advection problem, 
and the power law which characterises the divergence is 
universal. It should be noted that our studies incorpo- 
rate the important time-dependent correlations of the 
driving field. Moreover we show analytically that in the 
limit of a stationary interface, the equilibrium state of 
the sliding particles is an SCS, and determine the form 
of the scaling function. Remarkably, we find that the 
calculated equilibrium scaling function describes very 
well the SCS for the strongly noncquilibrium advection 
problem as well. 

The evolution of the one-dimensional interface is de- 
scribed by the KPZ equation [llj 
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where h is the height field, C,h is a Gaussian white noise 
satisfying ({ h (x,t){ h (x' ,*')) = 2D h S(x - x')5(t - t'). If 
the m th particle is at position x m , its motion is governed 
by 



dt 



dh 
dx 
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where the white noise Cm(t) represents the ran- 
domising effect of temperature, and satisfies 
((m(t)Cm(t')) = 2nS(t - t'). Equation © is a strongly 
overdamped Langevin equation of a particle in a poten- 
tial h(x,t) that is also fluctuating, with a determining 
the speed of sliding. In the limit when h(x, t) = h{x) 
is static, a set of nonintcracting particles would reach 
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the equilibrium Boltzmann state with particle density 
~ exp(— f3h(x)) at late times where j3 = a/ft. On the 
other hand, when h(x,t) is time dependent, the system 
eventually settles into a strongly noncquilibrium steady 
state. The transformation v = —dh/dx maps Eq. 
(with A = 1) to the Burgers equation, which describes a 
compressible fluid with local velocity v. The transformed 
Eq. (J2J describes passive scalar particles advected by 
the Burgers fluid. The ratio a/A > corresponds to 
advection, the case of primary interest in this paper, 
while a/A < corresponds to anti-advection (particles 
moving against the flow). 

Rather than analysing the coupled Eqs. and (0 
directly, we study a lattice model which is expected to 
have similar behaviour at large length and time scales. 
The model consists of a flexible one-dimensional lattice 
in which particles reside on sites, while the links or 
bonds between successive lattice sites are also dynamical 
variables which denote local slopes of the surface. The 
total number of sites is L. Each link takes values +1 
(upward slope — > /) and —1 (downward slope — > \). The 
rules for surface evolution are : for advection, choose 
a site at random, and if it is on a local hill (— > /\), 
change the local hill to a local valley(^ \/); otherwise 
leave it unchanged. For anti-advection, if a site is 
at a local valley, change it to a local hill; otherwise 
leave it unchanged [lSj. We use periodic boundary 
conditions, implying no overall tilt of the surface. After 
every N a surface moves we perform N p particle updates 
according to the following rule : we choose a particle at 
random and move it one step downward with probability 
(1 + K)/2 or upward with probability (1 - K)/2. The 
parameter K ranges from 1 (particles totally following 
the surface slope) to (particles moving independently 
of the surface). We define the ratio N s /N p as ui. The 
continuum Eq. l[T]l. valid for oj = 1, gets modified for 
u> ^ 1. While the first two terms on the R.H.S. of Eq. QJ 
get rescaled by u, the noise term gets multiplied by \fuj 
18] . The limit lu — > corresponds to the adiabatic limit 
of the problem where particles move on a static surface 
and the steady state is the thermal equilibrium state. 
In our simulations, we update the surface and particles 
at independent sites, reflecting the independence of the 
noises Ch(x,t) and ( m (t) EH- 

If ui is nonzero, the RMS displacement of each particle 
is proportional to t 1 ^ where z = 3/2 is the dynamical ex- 
ponent of the KPZ surface At large enough times 
t ^> L z , the system settles into a nonequilibrium steady 
state and develops interesting correlations, to which we 
now turn. We begin by discussing results for advection, 
for u> = K = 1. Figure Q shows the two-point (un- 
connected, scaled) density-density correlation function 
G(r,L) = (n,ini +r )L where rij is the number of parti- 
cles at site i. The total number of particles N = J2 n i 1S 
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FIG. 1: Two point scaled density correlation G(r, L) function 
(advection) for K = 1, o» = 2 (top curve), 1 (middle), 1/2 
(bottom). The line is a plot of Eq. Q with P = 4 - The 
lattice sizes are L = 256 (□, +, •), 512 (■, x, A) and 1024 
(°, *, A). 

taken to be equal to L. The numerical results give strong 
evidence that the scaling form 

G(r,L)~L- 6 Y(r/L) (3) 

is valid for r/L > with 9 ~ 1/2. An unusual point is 
that the scaling function Y(y) has a power law divergence 
Y(y) ~ y~ v as y — * 0, with v close to 3/2. Using this 
result, we can find the probability p(r) of two particles 
being at a distance r from each other in the limit L — + oo, 
p( r ) = ^rG(r,L) ~ -^kj2- This is in agreement with an 
exact result for two second class particles in the asym- 
metric exclusion process Another related quantity 
is defined by first dividing the system into L/l bins of 
of size I; then a particle is chosen at random and the 
number N(l,L) of the particles lying in the same bin as 
this particle is measured |6|. Using Eq.Q, we find that 
N(l,L) ~ ciL(l — C2l~ lJ+1 ), in better agreement with the 
numerical results for N(l, L) than the 1-independent form 
of |3 . 

The central curve in Fig.[5]shows the scaled probability 
P(n, L) that any site has occupancy n. This quantity, for 
n > 0, shows the scaling form 

with 5 = 1. As y — > 0, the scaling function f(y) seems to 
behave as with 7 ~ 1.15. We will see below that the 
values of exponents and the functional forms of the scal- 
ing functions in Eqs. and Q agree surprisingly well 
with those for the equilibrium case. Note that Eq. i fi} 
leads to (ri 2 ) = G(0, L) ~ L which is verified by direct 
simulation. From simulations we also find that the num- 
ber of occupied sites N occ = (1 — P(0, L))L varies as LP 
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FIG. 2: Scaled probability distribution P(n, L) for co = 
1/2, 1, 2 (K = 1). The line is a fit to Eq. ® with f3 = 2.3. 
The lattice sizes are L= 512 (o, x, □), 1024 (■, +, *). 



with <fi ~ 0.23, though the effective exponent seems to de- 
crease systematically with increasing L. A scaling anal- 
ysis yields 8 = v — 9 and <fi = d(-f — 2) + 1 with the latter 
holding for 1 > 5 > and 7 > 1. We also monitored the 
time dependent correlation function in the steady state. 
We find that (see Fig. |3J) the autocorrelation function 
G(t,L) = (n,i(0)n,i(t)) l scales with system size as 



G(t,L) 



(5) 



with z = 3/2. The scaling function shows a power law 
behaviour Y(y) ~ y~v with ip ~ 2/3 as y — ► 0. 

These results lead to a simple picture of a typical SCS. 
The scaling of the probability distribution P(n, L) with 
n/L and the vanishing of the probability of finding an 
occupied site (= N occ /L) suggest that a large number of 
particles (often of the order of system size) accumulate 
on a few sites; the scaling of the two-point density- 
density correlation function with L implies that the 
particles are distributed over distances of the order of 
L, while the divergence of the scaling function indicates 
clustering of large-mass aggregates. 

Now let us turn to the limiting adiabatic case u> —* 
corresponding to an equilibrium system of particles 
at inverse temperature (3 oc ln((l + K)/(l — K)) dis- 
tributed on a disordered, stationary surface. Instead of 
considering the number density of many noninteracting 
particles, we may consider equivalently the probability 
density of a single particle moving on the surface and 
average over all surface configurations, as in the Sinai 
model [lif . For the KPZ equation in one dimension, 
the distribution of heights in the stationary state is 



described by Prob[{/i(r)}] oc exp 



~2D h I ( 



dh(r') 
dr' 
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Thus, any stationary configuration can be thought of as 
the trace of a random walker in space evolving via the 
equation, dh{r) / dr = £(r) where the white noise £(r) has 
zero mean and is delta correlated, (£,(r)£,{r')) — S(r — r'). 
This is exactly the surface considered in the Sinai model. 
The probability p(r) = n r /L of finding the particle at 
position r is given by p(r) = cxp[— (3h(r)]/Z with the 
partition function Z = L L exp[— f3h(r')\dr' . 



The correlation function G(r, L)/L 2 = (p(ro)p(r + r )) 
involves an average over the surface configurations sam- 
pled from the stationary measure mentioned above: 



L- 2 G{r,L) = 



,-P[h(ro)+h(r +r)] 

z 2 



(6) 



The right hand side of Eq. © was evaluated exactly by 
Comtet and Texier m the context of one dimensional 
disordered supersymmetric quantum mechanics, where 
the right hand side of Eq. © is the correlation func- 
tion in the ground state wave function. In the scaling 
limit, r — > 00, L — > 00 with the ratio y = r/L fixed, one 
finds G(r,L) ~ L~ l / 2 Y{r / L) where the scaling function 
is given by 



Y(y) = 
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Note that the point r = is not part of the scaling func- 
tion; we have G(0, L) w (3 2 L/12. This formalism can also 
be used to calculate the equilibrium probability density 
P(p,L) [l^. Our results indicate that P(p,L) can be 
written as the sum of two parts : 



1 



ln 2 (L) 



(3 2 L 
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The first part refers to vacant stretches, and to the fact 
that the number of occupied sites occupies a vanishing 
fraction ~ (lnL) 2 /L of the system. The scaling function 
g(y) in the second part is given by 



g{y) 



y 



-K (y). 



(9) 



where K (y) is the modified Bessel function which 
has the asymptotic behaviour [— ln(y/2) — 0.5772...] as 
y — > 0. The theta function incorporates a lower cutoff 
on the validity of the scaling form and c is a constant 
of O(l). Thus we see that in the ui — > limit, the 
equilibrium state is an SCS. 

Surprisingly, these equilibrium results reproduce quite 
well the scaling exponents and scaling functions for 
G(r) and P(n) for n > 1 obtained numerically for the 
nonequilibrium case u) = K = 1, as can be seen in 
Figs, n an d [21 though with different values of /3. The 
correlation function matches with (3 ~ 4 while (3 ~ 2.3 



4 



CD 10 



X " a 
+ +. □ 

**x 

"Equilibrium" 



Nonequilibrium 



X 4 



X 



X X 
X*.. °»f 



0.01 

t/L z 



FIG. 3: Scaled density-density autocorrelation function 
G(t, L) (advection) in the nonequilibrium and equilibrium 
(adiabatic limit) cases. The dashed line shows y ~ x~ 2 ' 3 . 
The lattice sizes are L=256 (□, x), 512 (*, +). 



describes the probability distribution of number data 
well. However, P(0, L) (and thus N occ ) does not agree 
closely in the two cases. The equilibrium case can also 
be used to shed light on the dynamical properties of the 
nonequilibrium steady state. We compared our results 
for G(t, L) with the density-density autocorrelation 
function in the adiabatic w — > limit. To find the latter, 
we simulated a surface with height field h(r, t) evolving 
according to KPZ dynamics, and evaluated the density 
using the equilibrium weight p{r,t) — exp[— f3h(r, t)]/Z. 
As shown in Fig. the results with (3 — 4 agree with the 
autocorrelation function in the nonequilibrium system, 
apart from a numerical factor. 



small values of r/L (for G(r,L)) and n/L (for P(n,L)). 
Nonuniversal features emerge at large values of the 
arguments r/L and n/L. For u> < 1, we find that P(n) is 
nonmonotonic and shows a peak at large n (Fig.EJ). Nu- 
merically, we find that P(n = N)/P(n = N - 1) = 1/w 
(for u> < 1), which can be argued for by considering the 
processes which lead to the formation of an N particle 
cluster from an N— 1 particle cluster and vice versa |l8j . 

We have also studied anti-advection, where particles 
move opposite to the direction of surface motion. We 
again find that the density-density correlation function 
follows the scaling form of Eq. but with 6 = 0. 
The scaling function Y(r/L) diverges as {r/L)~ u with 
v ~ 0.33. Thus the steady state for anti-advection is also 
an SCS; the correlation function depends strongly on the 
size, and does not follow a size independent power law 
as in 7J. It is also interesting to study passive sliders 
on an Edwards- Wilkinson surface [l9( instead of the 
KPZ surface. Our simulations on this problem show an 
SCS having the scaling form of Eq. © with 9 = 0, and 
v ~ 0.66. As for advection, the relation S = v — holds 
in both these cases. 

In conclusion, we remark on two open problems. First, 
it would be worthwhile to understand better the cor- 
respondence we have noted between the SCS in the 
nonequilibrium advection and the SCS in the equilib- 
rium, Sinai problem, and in particular why different 
quantities in the former involve different temperatures 
in the latter. Secondly, the existence and extent of clus- 
tering induced by fluctuating potentials in higher dimen- 
sions remain to be studied and established. 



Why do the equilibrium results describe the non- 
equilibrium steady state so well? A partial explanation 
may lie in the fact that the statistical properties of 
the valleys where the particles settle down in both the 
cases are similar; under advection dynamics, particles 
slide down and reach the valleys relatively quickly, 
allowing them to explore the terrain. So, at least in a 
certain range of parameter space, the surface fluctuation 
noise mimics the effect of temperature insofar as it 
causes redistribution of particles over the landscape. 
Intriguingly, the temperature depends on the property 
under consideration — a phenomenon that deserves 
further study. 

The close correspondence between the SCS obtained 
in the nonequilibrium case and the equilibrium case 
u> = suggests there may be some universal features as 
parameters are varied. We confirm this by simulations, 
varying u> in the regime 1/4 < uj < 4. We find that for 
large L scaling is valid, with scaling functions which have 
universal exponents characterising singular behaviour at 
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